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We study the dynamics of the contact between a pair of surfaces (with properties designed to 
mimic ruthenium) via molecular dynamics simulations. In particular, we study the contact between 
a ruthenium surface with a single nanoasperity and a flat ruthenium surface. The results of such 
simulations suggest that contact behavior is highly variable. The goal of this study is to investigate 
the source and degree of this variability. We flnd that during compression, the behavior of the 
contact force displacement curves is reproducible, while during contact separation, the behavior is 
highly variable. Examination of the contact surfaces suggest that two separation mechanism are in 
operation and give rise to this variability. One mechanism corresponds to the formation of a bridge 
between the two surfaces that plastically stretches as the surfaces are drawn apart and eventually 
separates in shear. This leads to a morphology after separation in which there are opposing asperities 
on the two surfaces. This plastic separation/bridge formation mechanism leads to a large work of 
separation. The other mechanism is a more brittle-like mode in which a crack propagates across the 
base of the asperity (slightly below the asperity/substrate junction) leading to most of the asperity 
on one surface or the other after separation and a slight depression facing this asperity on the 
opposing surface. This failure mode corresponds to a smaller work of separation, those in which a 
single mechanism operates. Furthermore, contacts made from materials that exhibit predominantly 
brittle-like behavior will tend to require lower work of separation than those made from ductile-like 
contact materials. 



I. INTRODUCTION 

Many micro/nanoelectromechanicali systems 
(MEMS/NEMS) are based upon mechanical con- 
tacts that are only a few micrometers large. MEMS 
electrical switches must maintain high conductivity and 
mechanical reliability operating at radio-frequencies 
during their lifetime which may extend to years. How- 
ever, MEMS switches have been hindered by a lack of 
reliability. Extensive surface damage leads to failure of 
some MEMs devices after only several million open/close 
cycles.^ The resulting surface damage has been studied 
experimentally with atomic force microscopes in order to 
understand the effect of adhesion, thermal dissipation, 
and contaminationi^^— 

Although a variety of approaches have been used to 
study ideally flat surfaces ^^"^^ the contact surfaces of 
MEMS are rough^^ with a high density of nanoscale as- 
perities. Probabilistic fi^ii^ and fractali^"— approaches 
have been used to describe the correct topography of 
metal surfaces. MEMS contacts are intrinsically mul- 
tiscale, involving atomic bonding, defect, and fracture 
nucleation at the subnanometer scale, plastic and elastic 
deformation at the scale of a single nanoasperity as well 
as elastic deformation at the level of the entire MEMS 
switch. While the macroscopic elastic deformation in 
the presence of adhesion can be successfully character- 
ized in terms of Johnson, Kendall and Roberts (JKR) 
theoryi- and its more recent extensions ,■18--''' the plastic 
deformation and adhesion at the level of nanoasperity is 
difficult to treat analytically yet can be effectively mod- 



eled using molecular dynamic simulations. In this sense, 
a multiscale approach in which the microscopic elastic 
deformation of the contact is modeled by classical finite 
element simulations while the nanoscale contacts at the 
single asperity level are replaced by inelastic springs, with 
parameters taken from molecular dynamic simulations, is 
a promising avenue of MEMS research. 

A fundamental difference between mechanical behav- 
iors of microscopic MEMS contacts and single, nanoscale 
asperity contacts is that the latter may be governed by 
random events such as defect nucleation and fracture, 
while in the former these random events are likely to be 
averaged out due to the large number of asperities form- 
ing a macroscopic contact. Nevertheless, under some con- 
ditions, randomness at the nanoscale may lead to catas- 
trophic failure of the switch such as its permanent stic- 
tion. Thus, it is important to characterize the repro- 
ducibility and variability of the mechanical behavior at 
the single asperity level. 

In this paper, we focus on the statistical properties 
of the formation and breaking of the ruthenium single 
asperity contacts. Ruthenium contacts have proven to 
be more reliable than gold contacts, routinely surviving 
millions of cycles without significant degradation. The 
behavior of Ru single asperity contacts has recently been 
studied by molecular dynamic simulations and compared 
to that of gold contacts^^. In that work, a new embedded 
atom potentia l ^' -' —' was developed which accurately re- 
produces several key structural, thermodynamic and me- 
chanical properties of ruthenium, including its hexagonal 
close packed lattice structure, elastic constants, work of 



2 



adhesion, and stacking fault energy. It has been shown 
that while the behavior of Ru and Au asperity during 
compression is qualitatively similar, their behavior dur- 
ing separation is distinctively different. 

Au nanoasperity contacts^^"'^^ demonstrate typical 
ductile behavior. Namely, the asperity after sticking to 
an opposite substrate during compression, forms a sym- 
metric bridge, which gradually elongates on separation, 
while the neck of this bridge is getting thinner until its di- 
ameter becomes comparable to a single atom size. At this 
moment the contact breaks. In this case plastic deforma- 
tion occurs continuously without formation of a crack. 
This behavior is typical of ductile materials. On the other 
hand, a Ru nanoasperity contact was reported^^ to show 
a behavior which is more typical for brittle materials like 
glasses. Namely, the contact separation is characterized 
by crack formation and a sharp drop-to-zero of the tensile 
force during the unloading stage. 

In the present study we find that the behavior of Ru 
nanoasperity contacts is intrinsically chaotic. Depending 
on a slight change of initial conditions due to thermal 
noise, Ru bridges may break in brittle-like manner with a 
sharp drop of the tensile force or in a ductile-like manner 
which resembles the behavior of Au contacts. 

The brittle and ductile failure mechanisms are macro- 
scopic phenomena;^ while brittle materials like glass 
shatter very quickly, ductile materials like metals, can 
be deformed continuously. The origin of the mechanism 
is atomic in nature'^'^. At the atomic level the crystalline 
structure and the corresponding available glide planes for 
the motion of dislocations are the most important factor 
in the distinction between brittle and ductile materials. 
A non crystalline material, like glass, is brittle because 
dislocation movement is not possible. A crystalline ma- 
terial always has a certain degree of plasticity and the 
competition between a brittle and ductile behavior de- 
pends on how the interatomic bonds close to the tip of 
a crack respond to the locally large interatomic forces. 
Therefore at a microscopic level a material is more duc- 
tile or brittle depending on its crystalline structure and 
interatomic forces. A good operative definition of brit- 
tle and ductile behavior can be given by the separation 
mode. In Ref.^^, it was hypothesized that Ru contacts 
are more brittle than Au contacts due to the fact that the 
HCP lattice of Ru has fewer slip systems than the FCC 
lattice of Au and also due to the fact that Ru stacking 
fault energy (which is activated in plastic deformations) 
is 16 times larger than that of Au, while the Ru surface 
energy (which is activated during fracture formation) is 
only 3 times larger than that of Au. Our present study 
which shows the variability of the behavior of Ru contacts 
on separation does not contradict this hypothesis. 



II. METHODS 

We perform the molecular dynamics simulations using 
the LAMMPS package.^'* Pressure and temperature were 



kept constant using the Nose-Hoover thermostalj^ and 
barostat25 with a time step of St = 0.0025 ps. 

In the contact simulations, we create two substrate 
slabs facing one another. The substrate surfaces are fiat 
and parallel to the x — y plane and consist of 18 (0001) 
atomic planes of Ru hexagonal close packed (hep) lattice. 
Upon the lower substrate we place a homoepitaxial asper- 
ity, as shown in Fig. [T] The atoms of the top two atomic 
layers of the upper substrate and the bottom two layers 
of the lower substrate are not updated according to New- 
ton's equation of motion (as normal in molecular dynam- 
ics) but are displaced in the z-direction at constant rates 
in order to bring the surfaces into/out of contact and to 
compress/separate the substrates in z-direction. The x- 
and ?/-coordinates of the atoms in these layers are scaled 
to maintain zero net stress in the x- and y-directions. 
The system is periodic in the x- and y-directions. These 
boundary conditions are imposed to simulate large, finite 
substrates. 

The simulation cell contains a total of 62150 atoms 
of ruthenium, whose dimensions are = 97 A, Ly = 
103.1, and Lz = 121. 7A . The asperity is constructed 
as a homoepitaxial cubic island with dimensions 32.48, 
37.50, and 27.86 A (2891 atoms). Prior to the contact 
simulation the system is annealed in four stages, each one 
taking lOOps to be completed. First, the temperature is 
increased linearly from T = 3K to T = 1650K, close to 
the melting temperature of the ruthenium model used in 
this work, and subsequently reduced linearly to T = 3K. 
In the third stage, the system is heated again until it 
achieves T — 600K. Finally, in the fourth stage, it is 
equilibrated at T = 300K. 

After the asperity annealing is complete, we randomly 
assign atomic velocities from a Maxwell distribution at 
temperature T = 300K. The randomization of veloc- 
ities is done in order to examine the variability of the 
contact simulation results. Next, we displace the upper 
substrate toward the lower substrate at 0.07 A/ ps, while 
holding the lower substrate fixed (i.e., we hold the two 
atomic layers at the bottom of the lower substrate fixed) . 
When the distance between the top of the upper sub- 
strate and the bottom of the lower substrate reaches a 
predetermined limit, the sign of the velocity of the up- 
per substrate is reversed. The simulation continues until 
the upper and lower substrates are completely separated. 
The z-component of the force on the upper substrate is 
calculated during the simulation in order to determine 
the force-displacement relation. 

The embedded atom method (EAM) potential used 
to describe the interaction between Ru atoms may be 
written in the forni^ 

N-l N N 

f/=E E V{n,)+Y,F{p,), (1) 

where the subscripts i and j indicate each of the N atoms 
in the system, is the distance between atoms i and 
J, V{rij) is a pairwise potential, F{pi) is the embed- 
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ding energy function, pi = ^{rij) models the electron 
density at the position of atom i, and is another 

pairwise potential. The functional forms and parameters 
for V{rij), F{pi), and ^{rij) can be found in Ref. — . 
The interatomic potential adopted for the Ru atoms in 
this work reproduces the elastic constants and cohesive 
energy of Ru, and also gives stacking fault energy and 
surface energies that are in reasonable agreement with 
the experiment (for a comparison between experiment 
and model for several quantities of interest see Ref. -^). 
Stacking fault energies play an important role in plastic 
deformation and surface energy determines the work of 
adhesion; both are key elements in contact formation and 
separation phenomena. Despite the melting temperature 
of this model is low (1792 K) compared to the experi- 
ments (2607 K), this is not a problem since all simula- 
tions were performed at temperatures less than a quarter 
of the melting temperature of the model. 



III. CONTACTS 

Figure [1] shows a typical evolution of the Ru system 
during a single nanoasperity contact simulations. In this 
figure, we show two complete simulations corresponding 
to two different realizations of the initial atomic veloci- 
ties (exactly the same initial contact geometry). We can 
see that, depending on the initial conditions, the plastic 
behavior of the contacts during separation can be either 
ductile-like or brittle-like. In both cases in Fig. [H the 
system was subjected to the same degree of compression, 
and the upper and lower panels correspond to the same 
stage in the loading/unloading cycle. Figure [5] shows the 
force in the z-direction as a function of the displacement 
of the upper substrate for the two cases seen in Fig[T] 
[points (a)-(j) in this figure refer to panels (a)-(j) in Fig. 

The upper panels in Fig. [l] (a)-(e), show a case in 
which a crack forms near the base of the asperity and 
propagates across the asperity (nearly horizontally) as 
the two substrates are pulled apart [see Fig. [I{e)]. We 
characterize this behavior as brittle-like. The lower pan- 
els in Fig. [1] (f)-(j), show a nominally identical compres- 
sion/separation simulation, but with different, randomly 
chosen, initial atomic velocities. In this case, the behav- 
ior can be described as ductile-like in the sense that there 
is significant asperity necking (thinning) prior to separa- 
tion. We can see the difference between the brittle-like 
and ductile-like scenarios by comparing panels (d) and 
(i), and (e) and (j) in Fig. [TJ We see that in (d), the 
lower substrate has almost no contact with the asperity 
anymore and the tensile force drops to zero (see Fig. ^ . 
In contrast, in (i) the asperity forms a bridge between 
the two substrates. Gradual elongation of this bridge 
corresponds to the slow decrease of the tensile force be- 
tween the two substrates (see Fig. [5]). In this case, the 
tensile force persists to much larger separations than in 
the upper panel. Panels (e) and (j) show very different 



final configurations for the two cases. While in (e) we 
observe a horizontal crack at the base of the asperity, 
panel (j) presents two opposing asperities, which we call 
"stalactite" and "stalagmite". Stalactite/stalagmite- like 
final structures are commonly seen in ductile-like mate- 
rials such as goldi^ 

In Fig. [21 zero displacement corresponds to the first 
contact between the top substrate and the asperity [see 
Fig. mb) and (g)]. At large separation, which corre- 
sponds to negative displacement [see Fig. [T](a) and (f)], 
the force is zero (see Fig. [2]). As the two substrates ap- 
proach each other, the force becomes negative because of 
the short range attractive interactions between the atoms 
on opposing substrate surfaces (this distance is defined 
by the interaction range of the interatomic potentials). 
This attraction elastically stretches the two materials and 
there is a jump-to-contact^"— indicated by the first neg- 
ative spike in the force-displacement curve (e.g., note in 
Fig. [2] that the force becomes negative at z = 0). As 
the displacement increases beyond this point, the system 
starts to compress. There is an approximately linear rise 
in the force with displacement, punctuated by a series 
of relatively sharp drops. The linear increase between 
the drops corresponds to elastic compression. The sharp 
drops correspond to defect generation, migration and/or 
annihilation events. When the sign of the substrate ve- 
locity changes [panels (c) and (h) in Fig. [T], the system 
begins to recover. The unloading is initially character- 
ized by a long linear elastic region. Eventually, the force 
reaches a minimum (i.e., a maximum tensile force), fol- 
lowing which the tensile force slowly decreases to zero 
over a long displacement range. In this region, the over- 
all force-displacement trend is also interrupted by sharp 
jumps. During the ductile-like run (f)-(j) the two sub- 
strates separate at larger (negative) displacement than 
that at which the initial contact occurred (i.e., zero dis- 
placement). On the other hand, in the brittle- like run 
(a-e), separation occurs at nearly the same displacement 
as where the original contact occurred. 

Figure [3] shows the force versus displacement for differ- 
ent maximum displacements Zm (i.e., the displacement at 
which the sign of the substrate velocity switches). The 
cases analyzed in this work are Zm = 2.23, 5.78, 9.28, 
12.6, 16.28, 18.03, 21.53, and 23.28 A. For each z^, five 
independent runs were carried out with different distri- 
butions of initial atomic velocities. As seen in this figure, 
the force displacement curves show very similar behav- 
ior on loading but exhibit significant differences from one 
run to the next at the same Zm- This can be seen from 
the different elongations of the tensile parts of the force- 
displacement curves in Fig.[3l Some of these curves show 
a long tail from the minimum to the point of separation 
(this is what characterizes a material as ductile-like here) . 
Other curves with the same Zm show an abrupt separa- 
tion (a sharp approach to zero force). For those cases we 
refer to the system as brittle-like, since the sharp drop in 
the tensile force to zero is related to an abrupt separation 
of the asperity from one of the substrates. 



FIG. 1: The atomic configuration of the system for different displacements of the upper substrate (see also the following figure) 
as the two surfaces are brought together (loading) and separated (unloading). Panels (a)-(e) and (f)-(j) show the evolution of 
the system for two different distributions of the initial atomic velocities (see the text for details). The maximum displacement 
of the upper substrate is Zm = 21.23A. Panels (a)-(e) show a case in which the asperity behaves in a brittle-like manner, while 
panels (f)-(j) suggest ductile-like behavior. 




FIG. 2: Force versus displacement for Zm = 21.53A. A posi- 
tive force corresponds to pushing the surfaces together; nega- 
tive to pulling them apart. The zero of the displacement is set 
at the separation between surfaces corresponding to first in- 
teraction as the surfaces are brought into contact. The contin- 
uous line corresponds to brittle-like behavior, as represented 
by panels (a)-(e) of Fig. [T] and the dashed line corresponds to 
ductile-like behavior, as represented by panels (f)-(j) of the 
same figure. The points labeled as (a)-(j) here, correspond 
to the configurations shown in Fig. [1] denoted by the same 
letters. 



IV. CONTACT STATISTICS 

In order to compare the different properties for the 
same asperity contact system, we assign each run with 
the same Zm a unique identifier 1 < i < 5. We focus on 
the statistical properties of the contact dynamics of the 
Ru system by analyzing simulation results for five inde- 
pendent runs for eight different Zm- First, we describe 
the Zm — 23. 28 A case in detail. This case provides a rep- 
resentative summary of the overall qualitative features of 
the behavior for all Zm studied in this work. 

Figure m shows our results for — 23.28A for each 
of five runs, each one having different initial atomic ve- 




z(A) 

FIG. 3: Force versus displacement of the upper substrate to- 
wards the lower one for different maximum displacements, 
Zm- The curves show little variation on loading. We show 
five independent runs for each of eight maximum displace- 
ments: z^ = 2.23, 5.78, 9.28, 12.6, 16.28, 18.03, 21.53, and 
23.28 A. Each run begins from the same atomic configuration 
but with atomic velocities chosen at random from the same 
(Maxwell) distribution. 



locities. This figure shows that during the loading stage, 
the curves are insensitive to the initial atomic conditions 
of the system. After the maximum load is reached and 
the sign of the velocity of the upper substrate reverses, 
all curves experience a near linear drop, reaching a local 
minimum (maximum tensile force at z w 17A). Beyond 
this point, however, the curves begin showing consider- 
able deviation from one another. The sharp drop of the 
tensile force represented by the bold solid line (i = 1) is 
characteristic of brittle-like systems. On the other hand 
the doubled dot-dashed line {i = 5), the tensile force 
shows a long tail indicating significant necking (plastic 
stretching), which is a signature of a ductile- like mate- 
rial. The other three curves (bold dashed, thin dashed, 
and dotted, i = 2 — 4) correspond to intermediate cases 
between a brittle-like and a ductile-like behavior. 
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FIG. 4: Five independent runs for the case in which Zm ~ 
23.28A. The double dotted-dashed line (run i — 5) has a long 
tail, characterizing a ductile-like behavior. In contrast, the 
bold solid line (run i = 1) sharply approaches the zero-force 
level, indicating a brittle-like crack, Other curves indicated 
by the bold dashed line (run i — 4), thin dashed line (run 
i = 2), and dotted line (run i = 3) correspond to various 
intermediate behavior. See the text for more details. 



It is important to characterize the surfaces after com- 
plete separation of the two substrates, in order to un- 
derstand the morphological consequences of brittle-like 
and ductile-like behavior, as indicated in the force- 
displacement curves. Changes in surface morphology 
may have important consequences for the behavior of 
high frequency MEMS/NEMS switches which are sub- 
jected to repetitive contact formation and separation. 
Figure [5] show topographic maps (surface heights a func- 
tion of X and y) of the two opposite surfaces in the most 
ductile-like and brittle-like cases seen in Fig. U] (i.e., bold 
dot-dashed i = 5 and solid i = 1 curves). In the ductile- 
like case, both surfaces [Figs. [5] (a) and (b)] have large 
asperities with peaks in different positions. This final 
shape suggests that the metallic bridge between the op- 
posite surfaces shears apart. In contrast, after brittle-like 
separation we see a large asperity on one substrate op- 
posing a shallow indentation on the other substrate [Figs. 
[5] (c) and (d)] . This final shape suggests formation of 
a crack near the contact between the asperity and one 
of the two surfaces. This crack may propagate through 
the pyramidal defects which form during the compres- 
sion below the base of the initial asperity or deep inside 
the opposite surface deformed by the approaching asper- 
ity. Such defects were observed in Ruthenium contacts 
in Ref.^^. Interestingly, crack propagation may lead ei- 
ther to the asperity being retained on the substrate from 
which it originally came or be transferred to the oppos- 
ing substrate surface, depending upon which side of the 
initial asperity the crack formed. If the crack goes below 
the base of the initial asperity as in Fig. ^d), the entire 
asperity is transferred to the opposite side. In contrast, 
in other cases, the crack was observed to pass near the 
intersection of the asperity and the opposing surface. In 
this case, little material is transferred between the op- 



posing surfaces. 

In order to quantify ductile-like and brittle-like be- 
havior, we calculate the work of separation Ws, which 
graphically corresponds to the area between the separa- 
tion curve (negative F) and the F = line. Mathemati- 
cally this can be written as 

f^f(F=Q) 

Ws= Fdz' (unloading), (2) 

Jzi{F=Q) 

where (unloading) indicates negative dz while Zi{F — 0) 
indicates the point on the force-displacement at which 
F becomes for the first time equal to zero during the 
unloading as indicated on Fig. [S]and Zf{F — 0) indicates 
the displacement at the end of simulation at when the 
surfaces are completely separated and hence F — 0. Note 
that Ws > 0, since during unloading F{z) < for any 
z between Zf{F — 0) and Zi{F = 0) and the lower limit 
of integration, Zi(F = 0), is lager than the upper limit 
of integration, Zf{F — 0). We also compute the work of 
compression or loading 

Wc = / Fdz'(loading) -f / Fdz' (unloading), 

Jo J Zm 

(3) 

where (loading) indicates positive dz. Since Zm > Zi{F = 
0) but F{z) > 0, the second integral in (3) is negative. 
The results of the work of separation and the work of 
compression versus the maximum displacement z^ are 
shown in Fig. [T] From this figure, we see that, in gen- 
eral, the amount of energy necessary to separate the two 
surfaces increases with the maximum displacement Zm- 
However, the work of separation shows large fluctua- 
tions, corresponding to the difference in energy required 
to separate the surfaces in a brittle-like (lower Ws) or 
ductile-like (higher Ws) manners. Comparing Ws for var- 
ious runs with Zm — 23. 28 A with the force displacement 
curves (Fig. 2]), we see that indeed the curve with the 
longest tail {i — 5) corresponds to the maximal Ws, while 
the curve with a rapid fall in the tensile force to zero 
(i = I) corresponds to the minimal Ws- Other curves 
correspond to intermediate values of Ws ■ Comparing Ws 
for other z„i with the corresponding force-displacement 
curves, we find that the largest values of Ws correspond 
to the longest tails in the force-displacement curve, while 
the smallest values of Ws correspond to the most abrupt 
decrease of the tensile force. Thus Ws is a good overall 
indicator of the ductile-like versus brittle-like behavior. 
This is consistent with the definition of ductile- like versus 
brittle-like behavior in bulk materials. 

The error bars in Fig. [7] represent the standard devia- 
tion of the W^ data over the five runs for each maximum 



displacement: a — 



'a/N)Eli {Ws,-Ws)']'^\ 

where Ws — i^/N) J2iLi is the average over N — 5 
independent measurements of W. As the degree of com- 
pression (zjn) increases, the magnitude of the error bars 
grow. On the other hand, the relative values of a (i.e.. 
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FIG. 5: Topographic maps of the opposing surfaces after complete separation for Zm ~ 23.28A. Panels (a) and (b) illustrate 
the case of ductile-like separation corresponding to the doubled dot-dashed i = 5 line in Fig. [J] Panels (c) and (d) illustrate the 
case of brittle-like separation corresponding to the bold solid i = 1 line in Fig. U Panels (a) and (c) represent the final shape 
of the upper surfaces in Fig.[T]), while panels (b) and (d) represent the final shape of the lower surface (where the asperity is at 
the beginning of the simulation). The color bars indicate the height of the surface (in A), measured from the average surface 
height. Red and yellow regions correspond to parts of the surfaces protruding away from the bulk of the material, while dark 
blue colors represents depressions in the surface toward the bulk of the material. 




FIG. 6: Schematic representation of a typical force versus dis- 
placement curve in the loading-unloading process. The area 
of the hachured region is numerically equals to the work of 
compression whereas the area between the F{z) curve and 
F — line stands for the work of separation. Zi is the first 
point in which the force becomes zero after the maximum 
compression (zm) and Zf is the point in which F — again. 
See the text and equations ([2]) and ((3)1 for details. 



a/Ws) are nearly independent of z^, as seen from Fig. 
mja). Figure ISl^b) shows an analogous graph for the work 
of compression, which indicates that the relative fluctu- 
ation of the work of compression is only approximately a 



quarter of that for the relative fluctuation of the work of 
separation. Linear regression yields negligible regression 
coefficients in the work of separation and compression 
cases, 1.16 X 10~^A~^ and — 6.32x 10~^A~^, respectively 
(see the lines in Fig. [8|). 

We also analyze the number of asperity atoms which 
are transferred from the lower to upper substrate upon 
separation - see Fig. [9] For small degrees of compression, 
only a small portion of the asperity is transferred to the 
upper substrate. This is because at very low compres- 
sion, the deformation of the system is nearly elastic and 
the contact area is small. As the degree of compression 
increases, more asperity atoms become in contact with 
the upper substrate, increasing the degree of adhesion 
and contact area between the two surfaces - resulting in 
larger material transfer. We also observe that the fluc- 
tuation in the number of asperity atoms transferred to 
the upper substrate varies greatly, especially at large Zm 
(> 12. 6A), going from 10% to 90% of the asperity trans- 
ferred in some cases. 

Stalactites and stalagmites are formed on the oppos- 
ing surfaces during the separation process. The heights 
of these structures vary dramatically from run to run (at 
the same Zm)- One way to characterize the topography 
of the surfaces after separation is to compute the cross- 
correlation between the heights at the same (x, y) coordi- 
nates of these landscapes. In order to do this, we divided 
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FIG. 7: Work of (a) separation and (b) compression versus maximum displacement, Zm- The work of separation is the work 
expended to separate the two surfaces from one another, while the work of compression is the energy expended in compressing 
the system and then removing the applied force. The different symbols in (a) denote five simulation runs with different initial 
velocity distributions. In order to compare different properties of the same run, we assign each run (with the same Zm) a unique 
number i from one to five. The symbols o, □, o, A, and * correspond to i = 1,2,3,4,5, respectively. Ws (Wc) is numerically 
equal to the area between F{z) and the line _F = in the unloading (loading) part of the force-displacement curve. The error 
bars represent the standard deviation of the data shown for each Zm- 




FIG. 8: The standard deviation, ct, of the data in Fig. [7] divided by the work of (a) separation and (b) compression, versus 
2m. The quantity a /W is the relative deviation of the work over five different runs for the same Zm- This data shows that o is 
proportional to Ws and Wc for different Zm- The straight lines in the figures represent linear regression results. In both cases, 
the slopes are very small; 1.16 x lO'^^A"^ and —6.32 x 10~''A~^ for panels (a) and (b), respectively. The data for Zm = 2.23A 
were not plotted since the total work done was very small (nearly elastic deformation). 



the system into a 25 x 25 grid in the xy plane, in which 
each grid cell is a rectangle of dimensions Ax = L^jlh 
and Ay = Ly/2b. For each cell we find the atom with 
the minimal z — ~hu{x,y) on the upper substrate and 
the atom with maximal z — hi{x,y) in the lower sub- 
strate [where the surface profiles are h{x,y)]. Then the 
cross-correlation is computed as 



^ ^ hu{x,y)hi{x,y) 



Nn 



{x,yen} 



^ K{x,y) ^ hi{x,y), 



" {x..ym} 



where Vl is the set of cells which belong to the contact 
area at the point of maximal compression, and Nq, is 
the number of such cells. A^n typically goes from 136 
to 225, depending on the magnitude of Zm- The area 



iVf2 AxAy has the physical meaning of maximum contact 
area. We implement this condition in order to compute 
correlations only for the area directly affected by the con- 
tact. Our cross-correlation definition is devised in a such 
a way that the typical shape observed in the ductile-like 
samples [Fig. [3a),(b)], in which two asperities face each 
other, yields positive correlation, while the final shape re- 
sulting from brittle-like behavior [Fig.[SJc),(d)], in which 
an asperity on one side faces a depression on the other, 
yields negative correlation. Comparing Fig.fTOland Fig. [7] 
we see that the runs with the largest positive correla- 
tions always coincide with the runs with the largest work 
of separation. This observation is consistent with the ex- 
pectation that ductile- like behavior is associated with the 
formation of metallic bridges (between the opposing sur- 
faces) which elongate as the substrates are pulled apart. 
When the bridge finally breaks, the stalactite and the 
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FIG. 9: Normalized number of asperity atoms which remain 
attached to the upper substrate after separation (the number 
of atoms transferred to the upper substrate divided by the 
total number of atoms in the asperity, A'' = 2981) versus Zm- 
See the text for more details. 
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FIG. 10: Correlation between the stalagmites and stalactites 
heights on opposing surfaces versus the maximum displace- 
ment in the z-direction. See the text for details. 



stalagmite formed on the opposite surfaces facing one 
another. Brittle behavior is associated with crack for- 
mation which commonly propagates through the defects 
formed below the surface during compression, resulting 
in the asperity on one side opposing a depression on the 
other. 



DISCUSSION 



Effect of Initial Conditions 



500- 



-500- 




FIG. 11: Effect of the assigning different random velocities 
to the atoms prior to the annealing of the asperities so that 
the resulting asperities have slightly different shape. We plot 
force over the asperity atoms against displacement in the 
z-direction. Each curve corresponds to a different asperity 
prepared with the standard annealing time. We do not see 
any qualitative differences with the previously studied case in 
which the simulations differ only by the atomic velocities at 
the start of the compression stage [Fig. [4] 



were calculated for different starting nanoasperities. The 
force displacement curves were found to be drastically 
different for different nanoasperities, with the maximal 
tensile force varying by more than factor of two. 

The changes in initial conditions due to thermal noise 
are much more subtle but as shown in the previous sec- 
tion they have a strong effect during the separation stage. 
In order to address the question of how the initial condi- 
tions affect the results in the loading-unloading stage in 
greater detail, we perform five additional simulations in 
which we assign different random velocities to the atoms 
before the annealing of the asperity resulting in five dif- 
ferent initial asperities. We than assign to the atoms in 
each asperity the same random velocities and perform 
the loading-unloading cycle. 

We also perform five additional runs starting from five 
different asperities annealed for 200ps instead of the stan- 
dard annealing time of lOOps as describe in section HIl A 
longer annealing time guaranties an even larger struc- 
tural variability. Each asperity is annealed starting from 
different random velocities, while the same random ve- 
locities were assigned to all asperities before the loading- 
unloading cycle. 



The main question we address in this article is how 
the thermal noise alters the behavior of nanoasperity 
contacts. Clearly, large structural changes in the ini- 
tial asperity can lead to much greater variability than 
thermal noise. An example of such a study was given 
in Ref.^^ where the force displacement curves for multi- 
ple loading-unloading cycles of Au nanoasperity contacts 



The results for the standard annealing can be seen in 
Fig. [TT] while the result for the asperities with longer 
annealing time are shown in Fig. [TD In both cases, it 
is clear that the tail of the tensile force varies in a way 
similar to the results shown in Fig. [3l As expected the 
variability is more evident for the asperities annealed for 
a longer time. 



9 



500- 
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500- 



FIG. 12: Same as figure [TT] but for the two times longer an- 
nealing procedure. Force over the asperity atoms against dis- 
placement in the z-direction. Each curve corresponds to a 
diflerent asperity prepared with the longer annealing time of 
200ps. 
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x(A) 

FIG. 13: Comparison between the large system and the small 
one considered in section [III 



B. Effect of Asperity Size 



We investigate the effect of asperity size on tlie loading- 
unloading curves. We prepared a larger system consisting 
of 111664 atoms, with 7000 atoms in the asperity. The 
comparison between the structures of the large system 
and the smaller one considered in section Hill can be seen 
in Figure [T3l in which the upper substrate was removed 
for better viewing purposes. 

Five big asperities were prepared according to the stan- 
dard annealing procedure described in sectionlllland with 
five different initial random velocities. The same initial 
velocity distribution was used for the five asperities in 
the loading- unloading cycle. Here, only two maximal 
compressions were considered: 17.6 A and 25.0 A. The 
results are shown in Fig. [Til and Fig. [15] respectively. We 
find a high degree of variability, independently of system 
size and degree of compression. 



-500- 
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FIG. 14: Force over the asperity atoms against displacement 
in the z-direction for the large system. 
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FIG. 15: The same as in Fig. [TJlbut for maximal compression, 
25.0 A. 



C. Gold statistics 

Due to its popularity as a contact material, and as a 
complement to our previous work^^ we performed sim- 
ulation to understand the contact statistic of Au. We 
performed simulations of Au contact loading and unload- 
ing for different starting atomic velocities using the Au 
asperity generated as described in Ref.^^. Fig. [TOl shows 
the force-displacement curves for three different maxi- 
mum displacement. For each displacement we produced 
five runs assigning different initial velocities to the atoms 
before the loading-unloading stage. Like in Ru we do not 
find any variability in the loading phase. During unload- 
ing, a degree of variability is present, but most notably 
there is no change of the mechanism of separation. That 
is, in all cases we observe a ductile-like transition. To 
compare these results with those of Ru, we calculate the 
effective tensile strains Zs — Ws/\Fmt\ for five runs for 
Ru and Au, where Ws is defined in Eq. ^ and \Fmt \ is 
the maximal tensile force achieved in a specific run. Fig- 
ure [17] shows different values of for different maximal 
displacements for both materials. We see that the values 
of Zg for Au are typically three times larger than for Ru 
which emphasizes the much higher ductility of Au with 
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VI. CONCLUSION 




z(A) 



FIG. 16: Force versus displacement for Au contacts. Shown 
are 5 independent runs for three different maximum displace- 
ments Zmax- 




FIG. 17: Effective tensile strains Za — Wa/\Fmt \ for five runs 
and different maximal compressions 2max for Ru and Au, 
where \Fmt\ is the maximal tensile force achieved in a spe- 
cific run. 



respect to Ru. On the other hand the relative variability 
of Zs are almost the same for Ru and Au. For each Zm 
we find Cz, the standard deviation of Zg, and compare 
the relative standard deviation Ozj z^ for Au and Ru av- 
eraged over four different Zm- We see that this quantity 
for Au is 0.14 while for Ru it is 0.17. Thus the presence 
of two different modes of separation for Ru contacts does 
not make them more variable than Au contacts. 



We find substantial variability in the way two ruthe- 
nium surfaces separate after contact. This variability is 
a reflection of the different modes of separation available 
to the system (ductile-like versus brittle-like) and of the 
fact that these different modes are competing with one 
another. For the same degree of compression, simula- 
tions (that only differ in the initial thermal velocities of 
the atoms) can lead to either brittle-like or ductile-like 
behavior, as seen in Figs. [1][S1[71 andfTUl The ductile- like 
behavior can be described as the formation of a bridge be- 
tween the contacting substrates that narrows via plastic 
deformation and eventually separates via a shear process. 
This leads to the formation of asperities on both facing 
surfaces. The brittle-like behavior may be characterized 
as the propagation of a crack through the substrate in a 
region that is damaged by plastic deformation, resulting 
in an asperity on one surface facing a depression on the 
opposite surface. The work of separation is much greater 
for the case of ductile-like separation than for the case 
of brittle-like separation. While the work of separation 
(as well as its fiuctuations) gradually increases with the 
degree of compression Zm, the surface cross-correlation, 
as well as material transfer shows a dramatic increase 
of fluctuations for large compression (i.e., Zm > 12A 
or '--^ 3 atomic layers). Below this value, the material 
transfer is always small(constituting a small part of the 
initial asperity), while above this threshold the magni- 
tude of the material transfer becomes highly variable, 
ranging from 10% to 90% of the original asperity. Mor- 
phology cross-correlation analysis shows that in the large 
compression case, this variability can be attributed to 
the difference in separation behavior - ductile-like ver- 
sus brittle-like. The ductile-like contact separation mode 
seen in ruthenium is akin to the separation behavior seen 
in goldj^ On the other hand, the brittle-like separation 
mode observed here is uncommon in ductile materials. 
This makes ruthenium behavior very distinct from the 
behavior of gold. Further, contacts made from materials 
that exhibit predominantly brittle-like behavior will tend 
to require lower work of separation than those made from 
ductile-like contact materials. Some plastic deformation 
is, however, desirable, since good electrical contact per- 
formance requires large contact area and the availability 
of some plasticity will lead to an increase in contact area. 
In conclusion, although Ru contacts have two modes of 
separations, brittle-like and ductile-like they are more re- 
liable than Au contacts because the effective separation 
distance z^ for Ru contacts is several times smaller than 
for Au contacts and the variability of Zg and Ws for Ru 
and Au are almost equal. Large values of Zs for Au indi- 
cate formation of large stalactites and stalagmites on the 
opposite contact surfaces after separation, leading to the 
increasing roughness of the surfaces after multiple switch- 
ing cycles^^ and finally to large fluctuations of the stic- 
tion force and fast wearing of the contacts. In contrast, 
in Ru the asperities remain relatively small and although 
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some degree of uncertainty is present due to two possi- 
ble modes of separation, this uncertainty is comparable 
with that of Au contacts, which exhibit only ductile-like 
behavior. 
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